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(57) Abstract: Apparatus(20, 100) for assaying a target analyte in a localized tissue region (22) that may include the target and other 
analytes comprising: a light source (34, 104) that illuminates the region with light at each of a plurality of wavelengths at which 
light is absorbed and/or scattered by tissue in the region wherein light at a least one of the wavelengths is absorbed or scattered by 
the target analyte; a signal generator (40) that generates signals responsive to intensity of the light from the light source (34, 104) at 
different locations in the localized region (22); and a controller (32, 102) that: receives the generated signals; processes the signals 
to determine an extinction coefficient for light in the localized region at each wavelength; and determines the concentration of the 
target analyte responsive to a solution of a set of simultaneous equations having as unknown variables concentrations of a plurality 
of analytes in the region (22), one of which is the target analyte, wherein each equation in the set expresses a relatioship between the 
extinction coefficient, the absorption coefficient and/or the reduced scattering coefficient for light at a different one of the plurality of 
wavelengths and at least one of the equations expresses a relationship between the extinction coefficient and the reduced scattering 
coefficient. 
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PHOTOACOUSTIC ASSAY METHOD AND APPARATUS 

FIELD OF THE INVENTION 

The invention relates to non-invasive in-vivo methods and apparatus for determining 
the concentration of a substance in a body. 



Non-invasive methods for assaying a "target" analyte, such as for example glucose, 
comprised in a region of body tissue are known in the art. In a near infrared spectroscopy 
(NIRS) method, light at a plurality of different wavelengths in a near infrared band of 
wavelengths is transmitted into a tissue region of the body to assay a target analyte in the 
1 0 tissue region. Light at at least one of the wavelengths, a "target wavelength" is absorbed or 
scattered by the target analyte. Intensity of light at the different wavelengths that is transmitted 



through the tissue region or scattered out of the tissue region is measured. The measured 
intensities are used to isolate and determine the contribution of the target analyte to an 
absorption or scattering coefficient of the tissue region at the target wavelength in the 
15 presence of contributions to the absorption or scattering coefficient by other "interfering" 
analytes in the tissue. Known values for the absorption or scattering cross sections of the 
target analyte at the target wavelength and the determined contribution of the target analyte to 
the absorption or scattering component are used to assay the target component in the tissue. 



20 tissue that are averages over relatively long optical path lengths through the tissue of light 
used to acquire the measurements. As a result, NIRS methods and technologies generally 
suffer from poor spatial resolution. In addition, NIRS signals tend to suffer from noise 
generated by scattering of light at tissue interfaces, such as the skin, and tissue 
inhomogeneities. NIRS methods tend therefore to exhibit relatively poor signal to noise ratios. 

25 An article by G. Yoon, et al. "Determination Of Glucose Concentration in a Scattering 

Medium Based on Selected Wavelengths by Use Of an Overtone Absorption Band", in 
APPLIED OPTICS 1 March 2002; Vol. 41, No 7 describes an NIRS method and device for 
assaying glucose in a tissue medium. The method describes criteria for choosing discrete 
wavelengths for light used in performing an NIRS assay so as to reduce influence of 

30 interfering analytes in the tissue on the results of the assay. Whereas the method is based on 
measuring NIRS absorption spectra of the tissue medium, the tissue medium is assumed to be 
scattering as well as absorbing. The article describes a device for assaying glucose having a 
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BACKGROUND OF THE INVENTION 



However, NIRS methods provide concentration measurements of a target analyte in 
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light source and a detector that are used to measure absorption spectra for the medium that 
have their relative positions optimized so that "measured spectra can be independent of 
medium scattering". 

For many medical procedures it is advantageous to accurately determine concentration 
5 of a target analyte for tissue regions that are relatively spatially localized. For example, in 
assaying glucose levels for a patient it is generally advantageous to measure glucose levels in 
blood. To acquire such measurements, the measurements should be spatially localized to a 
blood vessel or blood vessels so that the measurements are not "diluted", for example, by 
glucose levels in interstitial fluids. NIRS methods and devices, because of their relatively poor 
1 0 spatial resolution generally cannot provide such localized assays. 

Methods of measuring concentration of a target analyte in a tissue region using a time- 
resolved photoacoustic effect or optical coherence tomography (OCT) can provide 
measurements resolved to a relatively high spatial resolution. 

In a method using a time resolved photoacoustic effect, light at at least one wavelength 
1 5 for which light is absorbed or scattered by the target analyte is used to generate photoacoustic 
waves in the tissue region. Pressure produced by acoustic energy from the photoacoustic 
waves that arrives at a suitable acoustic transducer or transducers is used to assay the target 
analyte at locations in the region at which the photoacoustic waves are generated. Locations at 
which the photoacoustic waves are generated can be determined to within about 10 microns 
20 axially, along a direction of propagation of the waves, and to within about 200 microns 
laterally. As a result, an assay of the target analyte can be spatially localized to relatively small 
volumes having axial dimensions of about 10 microns and lateral dimensions of about 200 
microns. 

In OCT, light from a semi-coherent light source comprised in an interferometer is split 
25 into a reference light beam and a light beam that illuminates the tissue region. Light from the 
reference beam is reflected from a mirror to an "interference region" in the interferometer 
where it interferes with light scattered from the tissue region that reaches the interference 
region. An interference signal in the interference region is generated substantially only for 
reference and scattered light that reach the interference region after traveling substantially 
30 equal optical path lengths. As a result, the interference signal is generated substantially only 
for light scattered from material located in a small volume of the tissue region for which the 
optical path lengths of scattered light and reference light are substantially equal. The 
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amplitude of the interference signal is substantially proportional to a scattering coefficient for 
material in the small volume and is used to assay the target analyte in the small volume. 
Optical coherence tomography can provide axial spatial resolution of about a micron and 
lateral resolution of about 3 microns. Spatial resolution of an assay provided by OCT assaying 
5 is therefore on the order of a small number of microns. 

However, absorption and scattering cross sections of a target and interfering analytes 
in a tissue region contribute to photoacoustic or OCT signals used to assay the target analyte. 
Accuracy of the assay is generally compromised if contributions to the signals from scattering 
cross sections of the analytes are not assessed and distinguished from contributions to the 

10 signals from absorption cross sections of the analytes. Prior art has not provided methods for 
assaying a target analyte in a tissue region responsive to photoacoustic or OCT signals for 
which scattering cross section contributions to the signals are assessed and distinguished from 
absorption cross section contributions to the signals. 

SUMMARY OF THE INVENTION 

15 An aspect of some embodiments of the present invention relates to providing assay 

apparatus that uses the photoacoustic effect to assay a target analyte in a spatially localized 
tissue region and accounts for scattering of light used to generate the photoacoustic effect in 
determining the assay. 

An aspect of some embodiments of the present invention relates to providing assay 

20 apparatus that uses OCT signals to assay a target analyte in a spatially localized tissue region 
and accounts for scattering of light used to generate the OCT signals in determining the assay. 

An aspect of some embodiments of the present invention relates to providing a method 
for incorporating the effects of scattering of light in assaying a target analyte in a tissue region 
using the photoacoustic effect and/or OCT. 

25 An assay apparatus in accordance with an embodiment of the present invention 

comprises at least one light source that illuminates the tissue region with light at each of a 
plurality of different wavelengths, hereinafter referred to as "mensuration wavelengths". For at 
least one of the mensuration wavelengths light is absorbed and/or scattered, optionally 
strongly, by the target analyte. 

30 In some embodiments of the invention, the assay apparatus comprises at least one 

acoustic transducer, which senses pressure in photoacoustic waves generated at different 
locations in the tissue region responsive to the light. Alternatively or additionally, the assay 
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apparatus comprises an OCT interferometer. Interference signals are generated by the 
interferometer between light scattered from material at different locations in the tissue region 
and a reference beam of light provided by the light source. An extinction coefficient for each 
mensuration wavelength is determined for the tissue region either responsive to signals 
5 generated by the at least one acoustic transducer or interference signals generated by the 
interferometer. 

The mensuration wavelengths are determined so that each extinction coefficient is 
dependent on concentration of at least one of a same plurality of "mensuration" analytes. One 
of the mensuration analytes is the target analyte and the remaining mensuration analytes are 

10 "interfering" analytes. Each extinction coefficient therefore defines an equation having as an 
unknown variable a concentration of at least one of the mensuration analytes. Together, the 
extinction coefficients define a plurality of simultaneous equations having as unknown 
variables concentrations of the mensuration analytes in the tissue region. The number of the 
plurality of mensuration wavelengths and therefore the number of simultaneous equations is 

1 5 equal to or greater than the number of the plurality of mensuration analytes. 

In accordance with an embodiment of the present invention, a wavelength dependent 
function, hereinafter a "scattering coefficient function", which is parameterized by at least one 
characteristic parameter, is used to provide a value for the scattering coefficient for at least 
one of the mensuration wavelengths that contributes to the extinction coefficient at the 

20 wavelength. The equation defined by the extinction coefficient for a given mensuration 
wavelength comprises a term which is the scattering coefficient function evaluated at the 
mensuration wavelength. An assay of the target analyte is provided responsive to a solution of 
the simultaneous equations. 

In some embodiments of the present invention, at least one characteristic parameter of 

25 the scattering coefficient function is determined from information extraneous to information 
used to determine the set of simultaneous equations. In some embodiments of the present 
invention, at least one characteristic parameter of the scattering coefficient function is 
determined from an extinction coefficient determined from signals provided by the at least 
one acoustic transducer or alternatively by the interferometer. 

30 In some embodiments of the present invention, the number of the plurality of 

mensuration wavelengths and therefore extinction coefficients is greater than the number of 
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the plurality of mensuration analytes and the simultaneous equations are used to determine at 
least one characteristic parameter of the scattering function. 

In some embodiments of the invention, the scattering coefficient function is 
determined assuming that optical scattering in the tissue region is Mie scattering. 
5 In some embodiments of the invention, the target analyte is glucose and an assay 

apparatus is used to provide in vivo measurements of glucose in a blood vessel in the body of 
a patient. 

There is therefore provided in accordance with an embodiment of the present 
invention apparatus for assaying a target analyte in a localized tissue region that may include 

10 the target and other analytes comprising: a light source that illuminates the region with light at 
each of a plurality of wavelengths at which light is absorbed and/or scattered by tissue in the 
region wherein light at at least one of the wavelengths is absorbed or scattered by the target 
analyte; a signal generator that generates signals responsive to intensity of the light from the 
light source at different locations in the localized region; and a controller that: receives the 

15 generated signals; processes the signals to determine an extinction coefficient for light in the 
localized region at each wavelength; and determines the concentration of the target analyte 
responsive to a solution of a set of simultaneous equations having as unknown variables 
concentrations of a plurality of analytes in the region, one of which is the target analyte, 
wherein each equation in the set expresses a relationship between the extinction coefficient, 

20 the absorption coefficient and/or the reduced scattering coefficient for light at a different one 
of the plurality of wavelengths and at least one of the equations expresses a relationship 
between the extinction coefficient and the reduced scattering coefficient. 

Optionally, the at least one equation that expresses a relationship between the 
extinction coefficient and the reduced scattering coefficient includes a dependence on the 

25 absorption coefficient. 

Additionally or alternatively the reduced scattering coefficient at at least one of the 
wavelengths is a measured value of the reduced scattering coefficient. 

In some embodiments of the present invention, the reduced scattering coefficient at at 
least one of the wavelengths is a value determined responsive to an analytic expression. 

30 In some embodiments of the present invention, the reduced scattering coefficient at at 

least one of the wavelengths is expressed as an analytic function. Optionally, the analytic 
expression is a function of at least one unknown variable having a value determinable 
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responsive to a solution of the simultaneous equations. Optionally, the at least one unknown 
variable is a concentration of at least one of the target analyte and the other analytes. 

In some embodiments of the present invention, the function comprises an expression 
of the form BAfC where X represents the wavelength and B and C are constants. 
5 In some embodiments of the present invention, the signal generator comprises at least 

one acoustic transducer that generates signals responsive to acoustic energy that reaches the 
transducer from photoacoustic waves generated in the region by the light. 

In some embodiments of the present invention, the signal generator comprises an 
optical coherence tomography device that receives light from the light source that is scattered 
10 from the region and generates an interference signal responsive to an interference pattern 
between the scattered light and light from the light source reflected by a reflector. 

In some embodiments of the present invention, the controller identifies and locates the 
localized region in a larger region comprising the localized region. 

Optionally, to identify and locate the localized region the controller: controls the light 
1 5 source to illuminate the larger region with light that is absorbed by a component characteristic 
of the localized region; receives signals generated by the signal generator responsive to 
intensity of the light from the light source in different locations in the larger region; uses the 
signals to assay the characteristic component in different localized regions in the larger region; 
and identifies and locates the localized region responsive to the assay. 
20 Optionally, the apparatus comprises at least one acoustic transducer controllable to 

transmit ultrasound, and to identify and locate the localized region the controller: controls the 
at least one transducer to transmit ultrasound into the larger region; receives signals generated 
by the at least one acoustic transducer responsive to acoustic energy reflected by features in 
the larger region from the transmitted ultrasound; and uses the signals to identify and locate 
25 the features and thereby the localized region. 

In some embodiments of the present invention, the localized region is a bolus of blood. 
There is further provided in accordance with the present invention a method of 
assaying a target analyte in a region of body tissue that may include the target and other 
analytes comprising: determining an extinction coefficient for light at each of a plurality of 
30 different wavelengths at which light is absorbed and/or scattered by tissue in the region and 
wherein light at at least one of the wavelengths is absorbed and/or scattered by the analyte; 
providing a value or an analytic expression for the reduced scattering coefficient at each 
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wavelength; and determining the concentration of the target analyte responsive to a solution of 
a set of simultaneous equations having as unknown variables concentrations of a plurality of 
analytes in the region, one of which is the target analyte, wherein each equation in the set 
expresses a relationship between the extinction coefficient, the absorption coefficient and/or 
5 the reduced scattering coefficient for light at a different one of the plurality of wavelengths 
and at least one of the equations expresses a relationship between the extinction coefficient 
and the reduced scattering coefficient. 

Optionally, the at least one equation that expresses a relationship between the 
extinction coefficient and the reduced scattering coefficient includes a dependence on the 
10 absorption coefficient. 

Additionally or alternatively determining the extinction coefficient at at least one of 
the wavelengths of the plurality of wavelengths optionally comprises: from a given location 
illuminating the region with light at the wavelength so as to generate photoacoustic waves in 
the region; determining a rate of decrease amplitude of the generated photoacoustic waves 
15 with increase of distance in the tissue region from the given location; and determining the 
extinction coefficient from the determined rate of decrease. 

In some embodiments of the present invention, determining the extinction coefficient 
at at least one of the wavelengths of the plurality of wavelengths comprises: from a given 
location illuminating the region with light at the wavelength; using optical coherence 
20 tomography to determine a rate of decrease of intensity of the light with increase of distance in 
the tissue region from the given location; and determining the extinction coefficient from the 
determined rate of decrease. 

In some embodiments of the present invention, the reduced scattering coefficient at at 
least one of the wavelengths is a measured value of the reduced scattering coefficient 
25 In some embodiments of the present invention, the reduced scattering coefficient at 

least one of the wavelengths is a value determined responsive to an analytic expression. 

In some embodiments of the present invention, the method comprises expressing the 
reduced scattering coefficient in at least one of the equations as an analytic function. 

In some embodiments of the present invention, the analytic expression is a function of 
30 at least one unknown variable having a value determinable responsive to a solution of the 
simultaneous equations. Optionally, the at least one unknown variable is a concentration of at 
least one of the target analyte and other analytes. 
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In some embodiments of the present invention, the analytic expression comprises an 
expression of the form BAfC where X represents the wavelength and B and C are constants. 

In some embodiments of the present invention, the method comprises identifying and 
locating the localized region in a larger region comprising the localized region. 
5 Optionally, identifying and locating the localized region comprises: illuminating the 

larger region with light that is absorbed by a component characteristic of the localized region; 
generating signals responsive to intensity of the light at different locations in the larger region; 
using the signals to assay the characteristic component in different localized regions in the 
larger region; and identifying and locating the localized region responsive to the assay. 
10 Additionally or alternatively, identifying and locating the localized region comprises: 

transmitting ultrasound into the larger region; generating signals responsive to acoustic energy 
reflected by features in the larger region from the transmitted ultrasound; and using the signals 
to identify and locate the features; using the identities and locations of the features to identify 
and locate the localized region. 
15 In some embodiments of the present invention, the localized region is a bolus of blood. 

BRIEF DESCRIPTION OF FIGURES 
Non-limiting examples of embodiments of the present invention are described below 
with reference to figures attached hereto, which are listed following this paragraph. In the 
figures, identical structures, elements or parts that appear in more than one figure are generally 
20 labeled with a same numeral in all the figures in which they appear. Dimensions of 
components and features shown in the figures are chosen for convenience and clarity of 
presentation and are not necessarily shown to scale. 

Fig. 1 schematically shows an assay apparatus assaying glucose using the 
photoacoustic effect, in accordance with an embodiment of the present invention; and 
25 Fig. 2 schematically shows an assay apparatus that uses both the photoacoustic effect 

and OCT to assay glucose, in accordance with an embodiment of the present invention. 
DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
Fig. 1 schematically shows an assay apparatus 20, hereinafter referred to as a 
"glucometer", assaying glucose in a "target region" 22 of a body part 24 of a patient, in 
30 accordance with an embodiment of the invention. Target region 22 is optionally located in a 
region 26 of soft tissue of body part 24 and comprises a body fluid, such as for example 
interstitial fluid, having a concentration of glucose. Optionally, target region 22 is a volume of 
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body fluid having a concentration of glucose and region 26 is a region of a fluid cavity 
containing the body fluid. For example, as in Fig. 1, target region 22 is a bolus of blood and 
the fluid cavity a blood vessel 23. 

Glucometer 20 optionally comprises a controller 32, a light source 34, optionally 
5 located in the controller, and an optic fiber 36 coupled to the light source. An end 38 of fiber 
36 is optionally mounted to a support structure 40, hereinafter a "probe head", to which an 
acoustic transducer or array of transducers is mounted Any of various appropriate acoustic 
transducers or array of transducers may be used in the practice of the invention. By way of 
example, in Fig. 1 probe head 40 has an array of acoustic transducers 42 positioned 
10 circumferentially around end 38 of optic fiber 36. Only two transducers of the array are 
shown. Probe head 40 is pressed to skin 44 of body part 24 to position end 38 of fiber 36 
close to or contiguous with the body part and to acoustically couple acoustic transducers 42 to 
the body part. 

To assay glucose in blood bolus 22 controller 32, optionally first, controls glucometer 

15 20 to locate blood vessel 23 and the bolus using any of various methods known in the art, such 
as methods described in PCT publication WO 02/15776, the disclosure of which is 
incorporated herein by reference. For example, to locate blood bolus 22 controller 32 may 
control transducers 42 to radiate ultrasound into region 26. Controller 32 processes signals 
generated by transducers 42 responsive to reflections of the radiated ultrasound from 

20 structures in region 26 to determine location of blood vessel 23. Alternatively, controller 32 
may control light source 34 to illuminate tissue region 26 with light that is relatively strongly 
absorbed by blood. Since the light is strongly absorbed by blood, photoacoustic waves are 
preferentially generated in blood vessel 23. Controller 32 processes signals generated by 
transducers 42 responsive to acoustic energy from the photoacoustic waves to image features 

25 in region 26 and locate blood vessel 23. 

In accordance with an embodiment of the invention, controller 32 then controls light 
source 34 to illuminate region 26 with at least one light pulse, represented by wavy arrows 50, 
at each of a plurality of mensuration wavelengths A^. The index i indicates a particular one 
of the mensuration wavelengths and satisfies the condition 1 < i< N\. 

30 The at least one pulse of light 50 at target wavelength X[ stimulates photoacoustic 

waves, schematically represented by starbursts 52, in tissue region 26 and bolus 22. 
Transducers 42 generate signals responsive to pressure in acoustic energy from photoacoustic 
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waves 52 that reach the transducers. The signals are transmitted to controller 32, which 
processes the signals in accordance with an embodiment of the invention, as described below, 
to determine glucose concentration in target region 22. 

In some embodiments of the invention, the at least one pulse of light 50 transmitted at 
5 different mensuration wavelengths X{ is transmitted at different times to illuminate bolus 22. 
In some embodiments of the invention the at least one pulse 50 comprises a train of pulses. In 
some embodiments of the invention the pulses in the train of light pulses at different 
mensuration wavelengths X[ are transmitted at different pulse repetition rates. Optionally, light 
pulse trains at different mensuration wavelengths are transmitted simultaneously. Signals 
10 generated by acoustic transducers 42 responsive to photoacoustic waves 52 that are stimulated 
by light pulse trains at different mensuration wavelengths are distinguished using signal 
processing techniques known in the art, such as appropriate heterodyning and phase locking 
techniques. 

Let intensity of a light pulse 50 transmitted at a mensuration wavelength X{ into body 
15 part 24, at a distance d from end 38 of fiber 36 be represented by l(k\,d). Assuming that d is 
larger than the mean free path for photons at wavelength \[, l(%i,d) may be written, 

(1) I(*i,d) = Io(^)exp(-cte(^i)d) 

where oe^i) is an extinction coefficient in the tissue of the body part for light at wavelength 
Aj and I^A^) is intensity of light in the light pulse at end 38 of fiber 36. The extinction 
20 coefficient is a function of an absorption coefficient o^^i) and a scattering coefficient cCs(Aj) 
in the tissue for light at the wavelength X{. Under the assumption of the diffusion 
approximation, the extinction coefficient may be written 

(2) Oe(fci) = [30^0(03^) + a's^j))] 1/2 
where, 

25 (3) = (l-g)a s (\i). 

a 's(^i) * s referred to as a reduced scattering coefficient and g is an anisotropy factor. 

Equation (2) may be rearranged to provide an expression for the absorption coefficient 

Oa(^i) 

(4) ota^) = 1/2 {-a' s (X0+[ a'siki^Hm)^^] 1/2 } . 

30 

At each wavelength X[, the absorption coefficient cc^A^) may be expressed as a sum of 
absorption coefficients of analytes in region 24 that absorb light at the wavelength X[. The 
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absorption coefficient of a given analyte is a product of an absorption cross-section of the 
analyte for light at wavelength X{ and concentration of the analyte in the body. Let the 
absorption cross-section of a "j-th" analyte at wavelength be represented by Oj(k[) and its 
concentration in blood bolus 22 by xj. In accordance with an embodiment of the present 
5 invention the mensuration wavelengths are chosen so that at each mensuration wavelength 
Ai, substantially only at least one of a same plurality of "N^" mensuration analytes contributes 
to Oa^i). One of the mensuration analytes Cj(X[) is glucose, the target analyte, and one of 
the Njl mensuration wavelengths is a target wavelength corresponding to the target analyte 
glucose for which target wavelength light is, optionally, strongly absorbed by glucose. The 
10 absorption cross section aj(Xj) for j=l and the wavelength X[ for i=l are arbitrarily assigned to 
represent respectively the absorption cross section for the target analyte glucose and the 
corresponding target wavelength. The absorption coefficient at wavelength X\ may therefore 
be written, 

N A 

(5) ot a (A d )= 2 a j(*>i>) = l/2{a' s fti)+[a' s (^^ 

j 

1 5 The Njl mensuration wavelengths provide a set of linear equations of the form of 

equation (5) in the unknown concentrations xj (1 < j < N^). The equations can be solved 
for any and all the mensuration analyte concentrations xj, and in particular for concentration 
x\ of glucose in blood bolus 22 if > and for each mensuration wavelength X{ of the 
wavelengths the extinction coefficient a e (A^) and reduced scattering coefficient a! s (Xi) are 

20 known. 

To determine the extinction coefficient a e (kj) for blood bolus 22 for each mensuration 
wavelengths A4 in accordance with an embodiment of the invention, controller 32 processes 
signals that transducers 42 generate responsive to pressure produced by photoacoustic waves 
52 at the transducers. 

25 Pressure sensed by acoustic sensors 42 responsive to photoacoustic waves 52 is time 

dependent. Pressure sensed at a time "t" following a time at which a light pulse 50 illuminates 
body part 24 arises from photoacoustic waves generated at locations in the body part for 
which distance "d" from acoustic sensors 42 is substantially equal to vt, where v is the speed 
of sound. (The transmission time of the light is negligible.) Let the pressure sensed by acoustic 

30 sensors 42 at time t responsive to a pulse of light 50 at wavelength X[ that illuminates tissue 
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region 26 be represented by P(^i,t). Then for photoacoustic waves generated at locations in 
tissue region 26 at a distance d from transducers 42, P (Jt^t) can be written: 

(6) P(^,t) = P(Xi,d/v) = KOa^ICM), 
where K is a constant of proportionality. 

5 Using equation 1, equation (6) may be rewritten, 

(7) P(?li,t) = P(Xi,d/v) = Ka^) {I o (Ad)exp(-0Ce(^i)d)} . 

From the time dependence of P(Ai,t), controller 32 determines which of the signals 
generated by transducers 42 responsive to P(hj,t) are generated responsive to photoacoustic 
waves originating at distances d from fiber end 38 corresponding to locations in bolus 22. 

10 (Distances d that correspond to bolus 22 are known from the location of blood vessel 23, 
which was determined as noted above.) From the signals responsive to photoacoustic waves 
52 originating inside bolus 22 controller 32 determines values for P(Aj,d) for a plurality of 
locations in blood bolus 22 at different distances d from end 38. The controller uses the 
determined values for P(Ai,d) and equation (8) to determine a value for oLq(Xi). Optionally the 

15 determined value for 0Lq(7^) is a best fit value that optimizes the fit of equation (7) to the 
determined values for P(^,d). 

To determine the reduced scattering coefficient o! s (k[) for mensuration wavelength X\, 
in some embodiments of the present invention, scattering of light in blood bolus 22 is 
measured at the wavelength. In some embodiments of the invention the scattering coefficient 

20 is determined from a wavelength dependent analytic function, i.e. a scattering coefficient 
function parameterized by at least one characteristic parameter. Optionally, the scattering 
coefficient function is determined assuming that scattering of light is substantially Mie 
scattering. As a result, as is known in the art, dependence of o! s (%i) on wavelength may be 
approximated by an expression of the form, 

25 (8) a' s (Xi) = EOfC 

In some embodiments of the present invention, values for the characteristic parameters 
B and C in equation (8) are determined for blood bolus 22 using methods known in the art, 
such as for example a method described in Mourant et al; "Mechanisms of Light Scattering 
from Biological Cells Relevant to Noninvasive Optical-Tissue Diagnostics"; Applied Optics 

30 Vol 37, issue 16, pg 3586-3593, June 1998. In addition, a reduced scattering coefficient 
a s(^R)» hereinafter a "reference scattering coefficient", for blood bolus 22 is determined for a 
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reference wavelength In terms of the reference wavelength and associated reference 
scattering coefficient, the reduced scattering coefficient o! s (X{) may be expressed by 

(9) a' s M = a' s ^R)(V^R)- C . 

In some embodiments of the present invention, a reference wavelength A,r for a tissue 
5 region is a wavelength for which the absorption coefficient oc^r) is known and the reference 
scattering coefficient o! s (kjt) is determined from equation (2) and measurements of an 
extinction coefficient (^(Ar) at the reference wavelength. In some embodiments of the 
present invention, the absorption coefficient CC^Ar) for a tissue region is known because the 
absorption coefficient of the tissue region is substantially determined by a component analyte 

10 whose concentration in the region is known. In some embodiments of the present invention, 
concentration of the component analyte is determined from a measurement of the extinction 
coefficient a e (Xj) for the region at a wavelength for which the extinction coefficient of the 
region is substantially equal to the absorption coefficient of the component analyte. 
Optionally, as described above, measurements of the extinction coefficient Ote(^R) are 

15 acquired from time dependence of signals generated by transducers 42 responsive to 
photoacoustic waves stimulated in the region. 

By way of a numerical example, for determining parameters of equation (9) required to 
determine a' s (Xj[) for blood, at 570 nm the magnitude of the reduced scattering coefficient 
oc' s (570) is between about 2 cm"* and about 3cm" 1. The magnitude the absorption coefficient 

20 ota(570) of blood at 570 nm is about 280cm" 1 . The extinction coefficient 0^(570) for blood at 
570 nm is therefore substantially equal to the absorption coefficient 0^(570) of blood. In 
addition, the absorption coefficient of blood at 570 nm is substantially equal to the absorption 
coefficient of hemoglobin. Furthermore 570 nm is an isobestic wavelength for hemoglobin at 
which the absorption cross-sections for oxygenated and deoxygenated hemoglobin are about 

25 equal. Therefore, at 570 nm the concentration of hemoglobin may be determined without 
having to know the ratio of oxygenated hemoglobin to total hemoglobin from a measurement 
of the photoacoustic effect at 570 nm. Equation (4) for blood at wavelength 570 nm becomes, 

(10) ota(570) = Oe(570) = c ah (570)x h , 

where 0^(570) is the absorption coefficient for hemoglobin at 570 nm and x^ is the 
30 concentration of hemoglobin in blood. In accordance with an embodiment of the present 
equation (10) provides a value for x^. 
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810 nm is another isobestic wavelength in the absorption spectrum of hemoglobin at 
which the absorption coefficient of blood is also dominated by the absorption coefficient 
0ah(81O) of hemoglobin. However, at 810 nm the reduced scattering coefficient is not 
negligible and equation (2) becomes, 
5 (11) Oe(810) = [3a ah (810)x h (a ah (810)x h + oc' s (810))]l/2 

Since x^ is known from equation (10), equation (11) may be solved to provide a value, 
in accordance with an embodiment of the present invention, for the reduced scattering 
coefficient a' s (810) at 810 nm. The scattering coefficient a! s (Xi) at wavelength X{ for blood 
may then be determined by using 810 nm for the reference wavelength A,r and oc f s (810) for the 
10 reference scattering coefficient in equation (9) to provide, 

(12) a f s (Xi) = a' s (810) (V810)" c . 

To determine the coefficient C in equation (12), optionally, the extinction coefficient 
Oq(X) is determined for at least two other, non-isobestic, wavelengths of light at which 
hemoglobin concentration substantially determines the absorption coefficient of blood. 
15 Suitable wavelengths are preferably wavelengths that straddle 810 ran, for example 950 nm 
and 700 nm. Let the straddling wavelengths be represented by X + and X", and collectively by 
Let the ratio of oxygenated hemoglobin to total hemoglobin in the blood be represented by 
S and the absorption cross sections for oxygenated and deoxygenated at wavelengths X- be 
represented by ^ahO^*) and a ahD(^) respectively, then the absorption coefficient cc^Ar 1 ) 
20 for hemoglobin in the blood at wavelengths A,* is, 

(13) ocah^i) = [OahQ^S + 0^(^(1 -S)]x h . 

Using equations (2), (12) and (13) and the extinction coefficients oLqQ^) determined 
for wavelengths S and the exponent C may then be determined from the two equations, 

(14) a e (X±) = [3aa h (?i±) ((Xah( X±) + tf s (810)(^±/810)-C)]l/2. 

25 Substituting the right side of equation (9) for a' s (^) in equation (5) provides an 

equation of the form, 

N A 

(15) ^oMx = l/2{a' s (X R )(X i /X R )-C+ WMHtyXjO-lC^4/3)a e (Ki)2]V2}. 

. J J 
J 

Equation (15) for the mensuration wavelengths 7^ provides a set of 
simultaneous equations in the unknown concentrations Xj, for each of which equations the 
30 right hand side the equation is known. In accordance with an embodiment of the invention, 
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controller 32 provides a value for the concentration x^ of glucose responsive to constraints on 
the concentrations xj defined by the simultaneous equations. Optionally, since water is a 
major component of living tissue and since the concentration of water is relatively labile at 
least one of the mensuration wavelengths is a wavelength, for example 1350nm, for which 
5 light is strongly absorbed by water and negligibly absorbed or scattered by other analytes in 
the body. 

For a number of mensuration wavelengths equal to a number of mensuration 
analytes, any of various well-known methods of manipulating and solving a set of 
simultaneous equations may be used to provide a value for x^. In some embodiments of the 

10 present invention a number of mensuration wavelengths is greater than N^, resulting in a 
number of simultaneous equations greater than the unknown concentrations xj. For such 
cases a suitable best-fit algorithm, such as a least squares algorithm may be used to provide a 
solution for concentrations xj. 

In some embodiments of the present invention, equation (15) is treated as an equation 

15 in (N^+2) unknowns, where in addition to the unknown concentrations of the 

mensuration analytes, the reference scattering coefficient a' s (A,R) and reference wavelength 
are considered to be unknown constants. Measurements of the extinction coefficient 
(Xe(Xi) are acquired for a plurality of mensuration wavelengths equal to or greater than 
(Na+2) to yield at least (N^+2) simultaneous equations of the form of equation (15). A set of 

20 at least (N^+2) simultaneous equation is sufficient to determine values for all concentrations 
xi, as well as for oc' s (A,r) and A,r. Controller 32 provides a value for the concentration xj of 
glucose responsive to constraints on the concentrations xj, reference scattering coefficient 
a 's(^R) anc * reference wavelength defined by the N% simultaneous equations. 

Similarly, in accordance with some embodiments of the present invention, the 

25 exponent "C n in equation (15), is also considered to be an unknown and measurements of 
oCe(Ai) are acquired for at least (N^ +3) mensuration wavelengths Xj. The at least (N^+3) 
extinction coefficient measurements yield at least (N^\+3) simultaneous equations of the form 
of equation (15). Controller 32 provides a value for the concentration x\ of glucose responsive 
to constraints on the concentrations xf, reference scattering coefficient cc' s (A,r), reference 

30 wavelength A,r and exponent C defined by the = (N^+3) simultaneous equations. 

In some embodiments of the invention for which the scattering coefficient is expressed 
as an analytic function having at least one unknown characteristic parameter which is a 
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concentration of at least one of the mensuration analytes. Optionally, the concentration of at 
least one of the mensuration analytes includes the concentration of the target analyte. If the 
analytic function representing the scattering coefficient at wavelength X is written S(A,,X), 
where X represents the set {xj} of concentrations of the mensuration analytes or a subset 
5 thereof, then equation (5) becomes, 

(16) oafli)- £a.(A.)x = ^{S^+tS^^^ote^jl^}. 

J 1 J 

J 

Similarly to the case of equation (5), the mensuration wavelengths provide a set of 
equations of the form of equation (16) in the unknown concentrations xj (1 < j < N^). 
The concentration of the target analyte is determined responsive to a solution of the set of 
10 equations. 

For some target analytes and conditions it is possible to choose an advantageous set of 
mensuration wavelength for determining concentration of an analyte in accordance with a set 
of equations of the form of equation (5) or equation 16. For example, as noted above at 570 
nm and 1350 nm the extinction coefficient for blood is substantially equal to the absorption 

15 coefficient of hemoglobin at 570 nm and water at 1350 nm respectively. At the isobestic 
wavelength 810 nm both the absorption coefficient and the reduced scattering coefficient 
contribute to the extinction coefficient for blood. It is possible and can be advantageous to 
assay glucose, in accordance with an embodiment of the invention, using these three 
wavelengths as mensuration wavelengths and hemoglobin, water and glucose as mensuration 

20 analytes. 

In particular, at 8 10 nm it can be advantageous to use a set of simultaneous equation at 
the mensuration wavelengths for which at least one equation has the form of equation (16) and 
the reduced scattering coefficient is represented by an analytic function S(^j,X). For example, 
if the concentrations of hemoglobin, water and glucose are represented by x^, x w and Xg, 

25 respectively, S(A^,X) in equation (16) optionally becomes S(Xj,xh,x w Xg), Optionally, 
S$i> x h> x w, x g) ma y b e expanded in a Taylor series to a desired order in the concentrations xfc, 
x w and Xg. Coefficients in the Taylor series may be determined from a suitable model and/or 
empirically. For example, the coefficients may be determined using an expression for the 
reduced scattering coefficient described in "Dynamic optical coherence tomography in studies 

30 of optical clearing, sedimentation, and aggregation of immersed blood"; Valery V. Tuchin, 
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Xiangqun Xu, and Ruikang K. Wang; APPLIED OPTICS Vol. 41, No. 1 , 258-271, January 
2002. Optionally, since at wavelengths 570 nm and 1350 nm the extinction coefficient is 
dominated by the absorption coefficient the reduced scattering coefficient is assumed to be 
zero and an expression for S(Aq,xj 1) x w Xg) is used only in an equation of the form (16) at 810 
5 nm. 

In the above examples, extinction coefficients (^(Aj) for the mensuration 
wavelengths and for the reference wavelength that are used to determine glucose 
concentration xj are described as being determined using the photoacoustic effect. In some 
embodiments of the present invention optical coherence tomography (OCT) is used to 
10 determine at least one of the extinction coefficients used to determine concentration of an 
analyte. 

OCT generally provides signals for determining an extinction coefficient having better 
SNR than photoacoustic effect signals at wavelengths for which the extinction coefficient is 
determined substantially by a reduced scattering coefficient. Photoacoustic effect signals 

15 generally have better SNR than OCT signals at wavelengths for which an extinction 
coefficient is dominated by an absorption coefficient. In accordance with an embodiment of 
the present invention, a glucometer for assaying glucose in a tissue region comprises at least 
one acoustic transducer and in addition an "OCT" interferometer. Photoacoustic signals 
generated by the at least one acoustic transducer are processed to determine extinction 

20 coefficients used to assay glucose for wavelengths at which an absorption coefficient 
dominates in determining a value for the extinction coefficient. Interference signals generated 
by the interferometer are processed to determine extinction coefficients used to assay glucose 
for wavelengths at which a scattering coefficient dominates in determining a value for the 
extinction coefficient. 

25 Fig. 2 schematically shows a glucometer 100, in accordance with an embodiment of 

the present invention comprising at least one acoustic transducer and an OCT interferometer. 
The components of the OCT interferometer are shown in a very schematic and simplified 
manner. Glucometer 1 00 is schematically shown assaying glucose in blood bolus 22 in blood 
vessel 23 of tissue region 26. 

30 Glucometer 100 optionally comprises a controller 102, and a light source 104, 

optionally located in the controller, that provides semi-coherent light at wavelengths for which 
it is desired to determine an extinction coefficient to assay glucose, in accordance with the 
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invention. An optic fiber 36 is coupled to light source 104 via an optical coupler 106. An end 
38 of fiber 36 is optionally mounted to a support structure 40 to which acoustic transducers 42 
are mounted. 

To assay glucose in bolus 22 controller 102 controls light source 104 to transmit at 
5 least one pulse of light into optical fiber 36 at each of wavelength for which it is desired to 
determine an extinction coefficient for bolus 22. An optical coupler 108 couples a portion of 
light transmitted by light source 104 along optic fiber 36 to an optical fiber 1 10 and transmits 
a portion of the light towards end 38 of fiber 36 from which end the light exits the fiber to 
illuminate tissue region 24. 

10 Some of the light that is transmitted along fiber 36 to exit the fiber at end 38 is 

absorbed in tissue region 26 and stimulates photoacoustic waves 52 in the region and some of 
the light is scattered by material in region 26. As in glucometer 20 acoustic transducers 42 
generate signals responsive to photoacoustic waves 52. Some of the scattered light reenters 
fiber 36 at end 38 and propagates back towards controller 102 through optical coupler 108. 

15 Light coupled to optical fiber 1 10 by coupler 108, exits the fiber from an end 1 12 and 

is reflected back into the fiber by a mirror 1 14. A portion of the light reflected back into optic 
fiber 1 10 is directed by optical coupler 108 to controller 102. When light reflected from mirror 
1 14 and scattered light from tissue region 26 that reenters optic fiber 36 reaches controller 
102, the light is directed by coupler 106 to a combiner 116. Combiner 116 superposes the 

20 scattered light and the reflected light at an interference region (not shown) to generate an 
interference signal. The position of mirror 114 relative to fiber end 112 is controlled by 
controller 102 to determine a desired path length from light source 104 to mirror 1 14 and back 
to combiner 116. The path length is determined so that substantially only light that is scattered 
in tissue region 26 from desired locations in the region generates an interference signal. Semi- 

25 coherent light source 104, couplers 106 and 108, optic fibers 36 and 110, mirror 114 and 
combine 1 1 6 cooperate to function as an OCT interferometer. 

Controller 102 processes signals generated by acoustic transducers 42 to determine 
which of the signals correspond to photoacoustic waves originating at different locations in 
bolus 22. Controller 102 controls the position of mirror 114 to scan bolus 22 and generate 

30 interference signals corresponding to light reflected from different locations in the bolus. 

For light at mensuration wavelengths having an extinction coefficient dominated by an 
absorption coefficient controller 102 optionally uses signals generated by transducers 42 
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corresponding to locations in bolus 22 to determine an extinction coefficient at the wavelength 
for the bolus. For light at mensuration wavelengths having an extinction coefficient dominated 
by a scattering coefficient, controller 102 optionally uses interference signals generated by 
combiner 1 16 corresponding to locations in bolus 22 to determine an extinction coefficient at 
5 the wavelength for the bolus. The controller uses extinction coefficients to assay glucose in 
the same way that glucometer 20 uses extinction coefficients to assay glucose. 

In the description and claims of the application, each of the verbs, "comprise" 
"include" and "have", and conjugates thereof, are used to indicate that the object or objects of 
the verb are not necessarily a complete listing of members, components, elements or parts of 

10 the subject or subjects of the verb. 

The present invention has been described using detailed descriptions of embodiments 
thereof that are provided by way of example and are not intended to limit the scope of the 
invention. The described embodiments comprise different features, not all of which are 
required in all embodiments of the invention. Some embodiments of the present invention 

15 utilize only some of the features or possible combinations of the features. Variations of 
embodiments of the present invention that are described and embodiments of the present 
invention comprising different combinations of features noted in the described embodiments 
will occur to persons of the art. The scope of the invention is limited only by the following 
claims. 
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CLAIMS 

1 . Apparatus for assaying a target analyte in a localized tissue region that may include the 
target and other analytes comprising: 

a light source that illuminates the region with light at each of a plurality of 
wavelengths at which light is absorbed and/or scattered by tissue in the region wherein light at 
at least one of the wavelengths is absorbed or scattered by the target analyte; 

a signal generator that generates signals responsive to intensity of the light from the 
light source at different locations in the localized region; and 

a controller that: 

receives the generated signals; 

processes the signals to determine an extinction coefficient for light in the 
localized region at each wavelength; and 

determines the concentration of the target analyte responsive to a solution of a 
set of simultaneous equations having as unknown variables concentrations of a 
plurality of analytes in the region, one of which is the target analyte, wherein each 
equation in the set expresses a relationship between the extinction coefficient, the 
absorption coefficient and/or the reduced scattering coefficient for light at a different 
one of the plurality of wavelengths and at least one of the equations expresses a 
relationship between the extinction coefficient and the reduced scattering coefficient. 

2. Apparatus according to claim 1 wherein the at least one equation that expresses a 
relationship between the extinction coefficient and the reduced scattering coefficient includes 
a dependence on the absorption coefficient. 

3. Apparatus according to claim 1 or claim 2 wherein the reduced scattering coefficient at 
at least one of the wavelengths is a measured value of the reduced scattering coefficient. 

4. Apparatus according to any of claims 1-3 wherein the reduced scattering coefficient at 
least one of the wavelengths is a value determined responsive to an analytic expression. 

5. Apparatus according to any of claims 1-4 wherein the reduced scattering coefficient at 

■ 

at least one of the wavelengths is expressed as an analytic function. 
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6. Apparatus according to claim 5 wherein the analytic expression is a function of at 
least one unknown variable having a value determinable responsive to a solution of the 
simultaneous equations. 

5 

7. Apparatus according to claim 6 wherein the at least one unknown variable is a 
concentration of at least one of the target analyte and the other analytes. 

8. Apparatus according to any of claims 4-7 wherein the function comprises an 
10 expression of the form BAr^ where X represents the wavelength and B and C are constants. 

9. Apparatus according to any of the preceding claims wherein the signal generator 
comprises at least one acoustic transducer that generates signals responsive to acoustic energy 
that reaches the transducer from photoacoustic waves generated in the region by the light. 

15 

10. Apparatus according to any of the preceding claims wherein the signal generator 
comprises an optical coherence tomography device that receives light from the light source 
that is scattered from the region and generates an interference signal responsive to an 
interference pattern between the scattered light and light from the light source reflected by a 

20 reflector. 

1 1 . Apparatus according to any of the preceding claims wherein the controller identifies 
and locates the localized region in a larger region comprising the localized region. 

25 12. Apparatus according to claim 1 1 wherein to identify and locate the localized region the 
controller: 

controls the light source to illuminate the larger region with light that is absorbed by a 
component characteristic of the localized region; 

receives signals generated by the signal generator responsive to intensity of the light 
30 from the light source in different locations in the larger region; 

uses the signals to assay the characteristic component in different localized regions in 
the larger region; and 
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identifies and locates the localized region responsive to the assay. 

13. Apparatus according to claim 11 or claim 12 wherein the apparatus comprises at least 
one acoustic transducer controllable to transmit ultrasound, and to identify and locate the 

5 localized region the controller: 

controls the at least one transducer to transmit ultrasound into the larger region; 
receives signals generated by the at least one acoustic transducer responsive to 
acoustic energy reflected by features in the larger region from the transmitted ultrasound; and 
uses the signals to identify and locate the features and thereby the localized region. 

10 

14. Apparatus according to any of the preceding claims wherein the localized region is a 
bolus of blood. 

15. A method of assaying a target analyte in a region of body tissue that may include the 
15 target and other analytes comprising: 

determining an extinction coefficient for light at each of a plurality of different 
wavelengths at which light is absorbed and/or scattered by tissue in the region and wherein 
light at at least one of the wavelengths is absorbed and/or scattered by the analyte; 

providing a value or an analytic expression for the reduced scattering coefficient at 
20 each wavelength; and 

determining the concentration of the target analyte responsive to a solution of a set of 
simultaneous equations having as unknown variables concentrations of a plurality of analytes 
in the region, one of which is the target analyte, wherein each equation in the set expresses a 
relationship between the extinction coefficient, the absorption coefficient and/or the reduced 
25 scattering coefficient for light at a different one of the plurality of wavelengths and at least one 
of the equations expresses a relationship between the extinction coefficient and the reduced 
scattering coefficient. 

16. A method according to claim 15 wherein the at least one equation that expresses a 
30 relationship between the extinction coefficient and the reduced scattering coefficient includes 

a dependence on the absorption coefficient. 
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17. A method according to claim 15 or claim 16 wherein determining the extinction 
coefficient at at least one of the wavelengths of the plurality of wavelengths comprises: 

from a given location illuminating the region with light at the wavelength so as to 
generate photoacoustic waves in the region; 
5 determining a rate of decrease amplitude of the generated photoacoustic waves with 

increase of distance in the tissue region from the given location; and 

determining the extinction coefficient from the determined rate of decrease. 

18. A method according to any of claims 15-17 wherein determining the extinction 
10 coefficient at at least one of the wavelengths of the plurality of wavelengths comprises: 

from a given location illuminating the region with light at the wavelength; 
using optical coherence tomography to determine a rate of decrease of intensity of the 
light with increase of distance in the tissue region from the given location; and 

determining the extinction coefficient from the determined rate of decrease. 

15 

19. A method according to any of claims 15-18 wherein the reduced scattering coefficient 
at at least one of the wavelengths is a measured value of the reduced scattering coefficient. 

20. A method according to any of claims 15-19 wherein the reduced scattering coefficient 
20 at least one of the wavelengths is a value determined responsive to an analytic expression. 

21. A method according to any of claims 15-20 wherein and comprising expressing the 
reduced scattering coefficient in at least one of the equations as an analytic function. 

25 22. A method according to claim 21 wherein the analytic expression is a function of at 
least one unknown variable having a value determinable responsive to a solution of the 
simultaneous equations. 

23. A method according to claim 22 wherein the at least one unknown variable is a 
30 concentration of at least one of the target analyte and other analytes. 
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24. A method according to any of claims 20-23 wherein the analytic expression comprises 
an expression of the form BX~^ where X represents the wavelength and B and C are constants. 

25. A method according to any of claims 15-24 and comprising identifying and locating 
the localized region in a larger region comprising the localized region. 

26. A method according to claim 25 wherein identifying and locating the localized region 
comprises: 

illuminating the larger region with light that is absorbed by a component characteristic 
of the localized region; 

generating signals responsive to intensity of the light at different locations in the larger 

region; 

using the signals to assay the characteristic component in different localized regions in 
the larger region; and 

identifying and locating the localized region responsive to the assay. 

27. A method according to claim 25 or claim 26 wherein identifying and locating the 
localized region comprises: 

transmitting ultrasound into the larger region; 

generating signals responsive to acoustic energy reflected by features in the larger 
region from the transmitted ultrasound; and 

using the signals to identify and locate the features; 

using the identities and locations of the features to identify and locate the localized 

region. 

28. A method according to any of the preceding claims wherein the localized region is a 
bolus of blood. 
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